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Abstract 

We consider a Randall-Sundrum two-brane cosmological model in the low en- 
ergy gradient expansion approximation by Kanno and Soda. It is a scalar-tensor 
theory with a specific coupling function and a specific potential. Upon introducing 
the FLRW metric and perfect fluid matter on both branes in the Jordan frame, 
the effective dynamical equation for the the A-brane (our Universe) scale factor 
decouples from the scalar field and B-brane matter leaving only a non-vanishing 
integration constant (the dark radiation term). We find exact solutions for the A- 
brane scale factor for four types of matter: cosmological constant, radiation, dust, 
and cosmological constant plus radiation. We perform a complementary analysis 
of the dynamics of the scalar field (radion) using phase space methods and exam- 
ine convergence towards the limit of general relativity. In particular, we find that 
radion stabilizes at a certain finite value for suitable negative matter densities on 
the B-brane. Observational constraints from Solar system experiments (PPN) and 
primordial nucleosynthesis (BBN) are also briefly discussed. 

Key words: two-brane cosmology, scalar-tensor theory in the Jordan frame, 
phase space methods, observational constraints 
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1 Introduction 



A great challenge for contemporary theoretical cosmology is to find and examine possible 
frameworks for explaining new observational data (dark energy, accelerated expansion, 
etc.) which may not be accommodated in the standard cosmological model based on the 
Friedmann solution of the Einstein equations. Investigations of alternative cosmological 
models include the braneworld scenario, where our visible Universe is one of the two 
hypersurfaces embedded in a 5-dimensional bulk spacetime, while the proper distance 
between the branes is measured by the radion field. 

The first phenomenological braneworld models were presented by Randall and Sun- 
drum |2] in order to solve the hierarchy problem of the Standard Model of parti- 
cle physics. A more fundamental theoretical setup for braneworlds arises in the 10- 
dimensional superstring theory (or the 11-dimensional M-theory) with 5 (6) compactified 
dimensions leaving us a 5-dimensional cosmological bulk spacetime [Sj. In models with 
several branes which move in the bulk, a collision of the branes can mimic the big bang 
[3] or the exit from inflation [S]. There are at least two approaches to get a 4-dimensional 
effective theory from the 5-dimensional theory. The first one, developed by Kanno and 
Soda [HI E|, is the gradient expansion method where the brane tension is taken to be 
much bigger than the energy density on the brane. The second approach, developed by 
Brax et al. [H], is known as the moduli space approximation. In this approach the brane 
positions are described by scalar moduli fields emerging from compactification assuming 
that the motion of the branes is slow. It can be shown that at least at the first order both 
approximations agree [0]. 

Following Kanno and Soda jjj, we consider the Randall-Sundrum type I cosmological 
scenario pQ with two branes (A and B) moving in a 5-dimensional bulk. Both branes 
are taken to be homogeneous and isotropic, and supporting energy-momentum tensors of 
a perfect barotropic fluid, with barotropic index T on the A-brane (identified with our 
visible Universe) and T B on the B-brane. The field equations of the effective 4-dimensional 
theory obtained by Kanno and Soda [7] are the equations of a scalar-tensor theory with 
one scalar field (interpreted as a radion) with a specific coupling function wfif). 

In recent decades, a lot of work has been done on studying scalar-tensor cosmologies 
with non-constant coupling function uu(^f). The general properties of the evolution of \1/ 
in the Einstein frame were investigated by Damour and Nordtvedt JU] and in the Jordan 
frame by Serna et al. [TT], while an analogous study in the case of a massive scalar field 
was presented by Santiago et al. [I2j. A comprehensive analysis of the early and late time 
cosmological models in the Jordan frame was carried out by Barrow and Parsons [T3] . 
Various observational constraints on scalar-tensor cosmologies were discussed by Santiago 
et al. Damour and Pichon ^5], Clifton et al. [16., and Coc et al. J7j, and for models 
with two branes by Palma [T%] . In these papers, most calculations are performed in the 
unphysical Einstein frame, although the Jordan frame is used for physical interpretations. 1 

1 For a lucid explanation of the relevance and relationship of the two frames, as well as a general review 



2 



The scalar-tensor cosmological model we consider in the Jordan frame stands out from 
those mentioned above in three respects. First, the second brane can contain matter and 
its matter tensor acts as an additional source in the equations of the gravitational and 
the scalar field. Second, the coupling function of the scalar field has a specific form 
u>(^) = 2(1-$) which is fixed by the construction of the model. Third, there are possible 
cosmological constants on both branes appearing via a specific potential of the scalar field, 
so in general we consider two-component source terms (matter + cosmological constants 
on both branes). 

The cosmology of this model was first investigated by Kanno et al. [2H] and later 
also by us [21] - 2 In this paper we obtain several new general analytic solutions for the 
flat (k = 0) A-brane scale factor, thus encompassing solutions for all interesting matter 
tensors: 1) cosmological constants and no matter, 2) radiative matter and no cosmological 
constants, 3) dust matter and no cosmological constants, as well as 4) radiative matter 
and cosmological constants. We also find the accompanying analytic solutions for the 
radion field in cases 2) and 4) above. We apply the methods of dynamical systems to 
study the evolution of the radion field, finding fixed points and drawing phase portraits 
related to the first three cases. The global picture emerging from this analysis is more clear 
and has some advantages over the usual procedure of studying the convergence of scalar- 
tensor theories towards general relativity, which relies on an analogy with a fictitious 
particle with velocity dependent effective mass [TUl HH H2] • Depending on the type and 
the (initial) density of A- and B-brane matter, the branes typically evolve towards infinite 
separation or to a collision, but for certain negative B-brane energy densities, the radion 
(corresponding to the inter-brane distance) is stabilized by an attractor at a certain finite 
value. Finally, we consider Solar system experiments and Big Bang nucleosynthesis to 
derive some observational limits on the parameters of the model. 

The plan of the paper is as follows. Sects. 2 and 3 give the basic equations of the 
model and clarify their general properties. Sect. 4 presents the analytic solutions, while 
Sect. 5 carries out a complementary analysis of the dynamics of the scalar field. Sect. 
6 follows with a brief discussion on observational constraints, and Sect. 7 contains some 
concluding remarks. 

2 Effective action and field equations in the Jordan 
frame 

The starting point is the 5-dimensional Randall-Sundrum action pQ 

J « s = ^|/ d 5 xV^(n -2A 5 )+ /dW=M£ m -£) + Jd A x^f(c B m -j: B ) , (i) 

of scalar-tensor cosmology, see the recent book by Faraoni • 

2 The cosmology and radion dynamics of RS two-brane model has been studied my many authors using 
various approaches [22| , the closest to ours being the moduli space approximation 1231 121] . 
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where 1Z is the scalar curvature, while h^ v and f^ v are the induced metrics on the A- and 
B-brane, respectively. We assume that the brane tensions can be divided into two parts 



S = (T + V 



<x n s + o B 



(2) 



with \a \ 3> | a |, | cr^ | ^> \a B \. The first part is balanced with the bulk cosmological 
constant A 5 , 



M 2 = K) 2 



6A, 



4 ' "°~ 

where i is the curvature radius of 5-dimensional bulk anti-de Sitter spacetime. The second 
part gives rise to unbalanced energy density on the branes which can be interpreted in 
terms of effective cosmological constants A, A s as 



6 

As = 



(3) 



Ac 



- A s c 2 

a 1 



° = ^r, ° B = ^-r- (4) 

Let us introduce a dimensionless scalar field called radion, which measures the 
proper orthogonal distance between the branes situated at y = and y = £ (at the fixed 
points of Si/Z 2 orbifold spacetime), 



d{x)= f ^dy = -Un(l-y), *6[0,1] 
Jo 2 



(5) 



in coordinates where = (x). Then, according to the low energy expansion scheme 
developed by Kanno and Soda, in a sense that the energy density of the matter on the 
brane is smaller than the brane tension |p/S| <C 1, |p s /£ s | <C 1, the 4-dimensional field 
equations on the A-brane can be written as [Zj 



G» v {h) 



K 



T\{h) + (l-^) 2 T%U) 



K 



St 



<r+ (1 - ^fa 1 



+1 fob _ s» + ^ffl ( - m\ a m 

^\ff V \ v v \<*J ^ \£2 \* 2 u 



2u + 3 
1 



T>) + (1-*)T*(/) 
do; 



4k 2 



2cj + 3 



[cx + (1 



2uo + 3 



I A* 



(6) 



(7) 



where ' denotes the covariant derivative with respect to the A-brane metric h^ v and k 2 = 
k§/£ is the 4-dimensional bare gravitational constant. The brane metrics are conformally 
related by /„„ = (1 — ^)h^ v while the function 



3 * 



(8) 



21 

describes the proper distance d between the branes as (cf. 

It is noteworthy that Eqs. (jHJ) and (0) coincide with the equations of the general 
scalar-tensor theory with a specific form (jHJ) of the coupling function uifi?) and an extra 



(9) 
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matter term from the B-brane. These equations can be derived from the 4-dimensional 
scalar-tensor action 



2k 2 



d xy—h 



VR(h) - ^p-¥ a ^ la - 2k 2 



2 r B 

m ' 



(10) 



where u(^>) is given by (jHJ) and V(^) denotes an effective potential introduced by cosmo- 
logical constants on the branes, 



V($) = a + a B (l - 



If C m = and a B = 0, the influence of the B-brane disappears from the action and we are 
left with scalar-tensor a theory with a constant potential V = a. The Eqs. © and (|7|) 
reduce to Einstein's general relativity when \I/ = 0, \? = 1, i.e., at u{^f) — > oo, ei — > oo. 



3 FLRW type two-brane cosmology 
3.1 Field equations 

Assuming that in the Jordan frame the Universe on the A-brane is described by the FLRW 
line element 

dr 2 



ds< 



-dt 2 + a 2 (t) 



+ r 2 (d9 2 + sin 2 0d<f) 



k = 0,±1 



;i2) 



and the matter on the branes is modeled by a perfect barotropic fluid with respect to the 
corresponding brane metric 



T. 



[IV 



(p + p)Uf,u u + ph 



p = {Y-\)p. 



T. 



(p B +p B )u B u B +p B f IMU , p b = (r - i)p s , 

we can write the 4-dimensional field equations (jHJ) and (J7|) as follows: 

^2 



K 2 V K 2 p K 2 (1 - *) 

* ■ 4*d - + y"^ + + y * p 



2 B _A 

a 2 



(13) 
(14) 

(15) 



2# + 3H 2 



^- ¥ (i-*)V + ¥ ^ 



* * 3 ^ 2 fc 

* " 2i/ ? " 4*(l-tt) ~ ^ ' (16) 



-3#* - - 

.2 



2 (1 - *) 3 V 



+y (1 - *) (p - 3p) + y (1 - tf) V - 3p^ 



(17) 



In the case \j? = 0, \& = 1 these reduce to the general relativistic FLRW equations with 
cosmological constant and barotropic perfect fluid. Note that since cosmological constant 
is included separately, its modeling by barotropic index T = is redundant. 
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The Bianchi identity yields the following equations (conservation laws) for the matter 
fluids: 



a N~ 3r 



p + 3HTp = p = Po[ — ) (1? 



- 3 r B 



™ r v-^(rV v=0 - " a ="Hv^^J (19) 

(here we assume that p/p > 0, p B / p B > 0, allowing the possibility of exotic negative 
energy density on the B-brane). Combining Eq. (jTHj) with Eq. (jTHj) we get a relation 
between the energy densities on both branes 




(20) 

The meaning of integration constants is given by Eqs. (fTHj) and (jl9j) : at a fixed moment 
t we have a(t ) = a , p(t ) = p , p B (t ) = p B , *(t„) = tf . 

Combining Eqs. (fT5|) and (fT6|) and using the radion equation (JT7J) we get the dynamical 
equation for if, 

H + 2H 2 = ^ 2 a+^(p-3p)-^, (21) 

which does not contain any additional terms compared with the usual FLRW cosmology. 
This happens due to the specific form of the coupling function (JHJ) and the potential (fTTf 
which cancel all terms containing radion, B-brane matter and B-brane cosmological con- 
stant. Yet, the differences from the general relativistic FLRW model are concealed in the 
bare gravitational constant n 2 which does not coincide with the Newtonian gravitational 
constant (see Sect. 6), and also in the form of the constraint equation (|15j) . 
In the case p ^ 3p (T ^ |) the first integral of Eq. (|2*Tj) reads 



d(a 2 



da 4 + po af a^ 3r -^ka? + Caf 



o 




(22) 



where C is a constant of integration. Using the conservation law (|18|) it can be written 
as an expression for the Hubble parameter, 



Kj k, / cl \ k Kj C 



H 2 = -a + - Po (- + _ . (23) 

3 3 \a J a z 3 \a J 

The last term is forbidden in the case of the standard FLRW cosmology due to the 
Friedmann constraint but is allowed in the present case. Let us note that Eq. (J2*3*j) 
coincides with the Friedmann equation of the fine-tuned one-brane cosmological model 
j2H] i n the low energy approximation where the term ~ p 2 is neglected. The integration 
constant C is geometrical in nature and can be interpreted as the value of the energy 
density of the dark radiation on the A-brane at the moment to- 
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In the case p — 3p (T — |) an equation corresponding to Eq. reads 



& = ^ + ^(-1 ■ (24) 



k 2 k k 2 K ( a 
3 a 2 3 \a , 

where K is an integration constant which can be redefined as a sum of the initial value 
of radiative matter density po an d the initial value of the dark radiation density C, K = 
Po + C. For k = 0, a = the exact analytic solutions for a(t) and are presented in 
Subsect. 4.2. 

Eqs. (|2"Tj) . (J2"3~j) imply an expression for acceleration, 



a k 



2 



((2-3r)p + 2a)-— ^, (25) 

which differs from the corresponding expression of the FLRW cosmology only due to 
the last term. To get an additional drive for acceleration besides the usual cosmological 
constant, we see the integration constant C must be negative, C = —\C\. 3 However, 
in an expanding Universe the last term diminishes quickly in time, so it could have been 
influential in the early Universe, but is negligible at present. 

3.2 Two-brane cosmology in the conformal time 

In order to clarify some general properties of the model let us introduce the conformal time 
on the A-brane, dt = a(t)dr] as usual, and let us redefine the scalar field as 1 — ^ = x 2 - 
Then from Eqs. ()15)) -([17 p the following equations follow: 

(a') 2 - ((a X )') 2 + k{a 2 - {a X f) = ^a\a + o* X " + P + pV) , (26) 

K 2 

a" + ka = — a 3 (4a+ (4 - 3I» , (27) 
6 

6~ 

with a prime here denoting d/drj. The result is not surprising taking into account that 
the B-brane metric is related to the A-brane metric by a conformal transformation with 
the conformal factor \ 2 . 

Dynamical equations (J27j) and (J2H|) can be integrated and their first integrals read 

(a') 2 = ^a 4 (a + P )-ka 2 + ^C, (29) 

{{a X )') 2 = -j(a X ) 4 (a° + p°) - k{a X f + C* , (30) 

where C B is another constant of integration. Eq. (|2l)|) constrains it to be C B = C. We 
see that the influence of the B-brane on the A-brane manifests itself only in permitting 



(a X )" + k(ax) = -Uax?{^ + {A-3T")p*) , (28) 



3 The same phenomenon was noticed before in the studies of a single brane |27| . 
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non-vanishing values of C which must be proportional to the value of the integration 
constant C B in the solution for scale factor (ax) of the B-brane. If C = the cosmologies 
on the both branes decouple and amount to the usual FLRW models. A non-vanishing C 
can be interpreted as an initial value of the energy density of the dark radiation on the 
A-brane (see Subsect. 2.2) 

4 Analytic solutions for a flat (k = 0) Universe 

In general, Eqs. (|15 |) - (|17|) describe a cosmological model with a two-component source 
term containing the cosmological constant and a perfect fluid matter (r 7^ 0) on both 
branes. Cosmo logically interesting matter tensors are radiation (r = 4/3) and dust 
(r = 1). However, the full analytic solutions of Eqs. ([15 |) -([17 )) can be found only in some 
special cases. 

4.1 Cosmological constant dominated Universe 

In the case of vanishing matter density (p = 0,p = 0) and positive cosmological constant 
on the A-brane (a > 0) the solution of Eq. (|21|) for a flat Universe (k = 0) reads 4 



th (±2 (t - tx)) \ a\ ch (±2 (t - t x 

H = I \ \ , ==> a 2 = { V 

±V^ cthf±2^/^ (t-h)) [ a\ sh (±2^/^ (t-i 

(31) 

where a%, t\ are integration constants. The solutions 1)31)) approach asymptotically in 
time the de Sitter solution if 2 = ycr, a > 0. In an expanding Universe (if > 0) the 
solution H ~ th(i) is accelerating (H > 0), while H ~ cth(i) is decelerating (H < 0). The 
influence of the dark radiation can be seen from the first integral ()23j) where integration 
constants are included differently, 

* = ,32) 

Here for the expanding Universe the last term on RHS is quickly decaying away leaving 
simply de Sitter expansion independent of dark radiation. 5 

For negative cosmological constant, a < 0, the corresponding solution reads 



' K 2 \a\ / K 2 \a\ 




cos 2\ ±— ^ (t-t 



(33) 

It describes a Universe which contracts into a singularity in finite time. 



4 The solution with th was previously written down in the moduli space approximation |24| . 
5 In this case (C = 0) an analytic solution also for the scalar field was presented in ref. |21| . 
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4.2 Radiation dominated Universe 

In the case k = 0, T = 4/3, o = Eq. (|21|) acquires the general relativistic FLRW form, 

H + 2H 2 = 0, => l__L = 2 (t-t 2 ), (34) 

where i?2, £2 are constants of integration, #(£2) = #2- The Hubble parameter decreases 
(H < 0) and the scale factor is 



a = a 2 ^J2H 2 (t - t 2 ) + 1 , (35) 

where a 2 is a constant of integration, a(t 2 ) = a, 2 . From Eqs. (J53j) and (|3*3|1 the Hubble 
parameter can be given in terms of the scale factor, 

tt2 „4 

H 2 = • (36) 

a 4 

Since it must coincide with the first integral (J24)) the initial value C of the dark radiation 
density at to can be expressed in terms of integration constants, 

C + Po = —^. (37) 

Kj CLq 

Eq. (|T7j) for a redefined radion field x 2 = 1 — ^ takes a simple form, 

X + 3#x = 0, (38) 

and its solution reads 

X - Xoo = ; 92 =— , (39) 



H 2 y/2H 2 (t -t 2 ) + l 

where Xoo and g 2 are integration constants. The constraint equation ()15j) with k = 0, 
V = allows us to determine one of them 

^oo = 1 - xl = ^4 (Po + pEXo) 4 , (40) 

where po> Pen Xo> a o are initial values at t = to m the solution of the conservation laws 
(|18|) . (|19p. In a realistic cosmological model, the radiation dominated era has a finite 
duration t € [tin, tout]- Substituting t 2 = t in and if 2 = H in into Eqs. (J35J) and (J59*j) . the 
change in the scalar field can be given in terms of e-folds, N = hx(a out /a in ), as 



_ Xin ( 1 „— iV 

Xout Xin 



(1 - e^) . (41) 



If the duration of the radiation dominated era is at least several e-folds, the change of the 
scalar field is approximately determined only by the values at tj n . This also means that 
after a few e-folds the scalar field remains approximately constant. 

The constraint between integration constants (|3U|) can also be used for putting Eq. 
f!36|) for H 2 into a form suitable for comparison with general relativity, 

# 2 = 7^ (l + (1 - *o) 2 ^ )p. (42) 
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Denoting the Hubble parameter of general relativity by H G and taking into account the 
standard Friedmann equation = (8nG/3)p (here G denotes Newtonian gravitational 
constant which in the context of scalar-tensor theory depends on the asymptotic scalar 
field, see Subsect. 6.1) the speed-up factor £ for radiation dominated Universe turns out 
to be constant, 

^5^( i+(i -*° )2 f)- (43) 

In terms of the initial values of radiative energy densities po, C, it takes a simple form 

£ 2 = ^(l + -V (44) 

Combining the last two equations, one may also express 

c = p (l-^ 00 )+pg(l-^o) 2 

This tells that having dark radiation contributing towards acceleration, i.e., C < 0, is 
only consistent with negative energy density for B-brane radiation, po < — (fzfffi Po- Yet, 
since the RHS of Eq. (|37jl is positive, it is clear that the net contribution of radiation and 
dark radiation to acceleration is always negative. 

4.3 Dust dominated Universe 

Integrating the first integral (J23j) with k = 0, a = 0, T = l leads to a cubic equation for 
the scale factor which has one real and two complex solutions. The real solution reads 



il,i C 2 ,i C \ 
«3 ( « + 2-j b-s + - ) , (46) 
A Po Po, 



where t 3 , a 3 are integration constants and the dimensionless quantity b can be written as 
b = p- 3 (z K 2 p* {t - hf - 8C 3 + y/9 K y (t - t 3 ) 4 - 48CWp$ (t - t 3 ) 2 ) • (47) 

The meaning of t% is clarified by the relation a(t 3 ) = 0*3^, which also implies C > 0. 
Note that in the case of vanishing dark radiation there is a singularity at the moment £3. 
If C 7^ and we take to to be the present time, then t 3 < t , if at present ^ < 1, and 
t 3 > t , if — > 1. 

The Hubble parameter as a function of the cosmic time reads 



f^Z 2(l-4(^) 2 6"l 
H=\ Po . V p0 '- . (46 



3^ P o {t - t 3 ) 2 - 16f (l + 2(£ 6-1) + 4 6-1 



In the case of vanishing dark radiation C = we get the FLRW value for the Hubble 
parameter H = |(t — ts)~ l . 
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4.4 Universe with cosmological constant and radiation 

Eq. (|21|) with = 0, a / does not depend on radiation matter with T = |, so 
its solutions coincide with those for cosmological constant dominated Universe ()31|). In 
particular, Eq. (f3*2*|) holds and the late time Universe is dominated by the cosmological 
constant o alone. 

For getting an equation for the radion field x 2 = 1 — ^ let us solve the corresponding 
constraint equation (fi3j) as an algebraic equation for x 



K 2 



X = -H X ±^H 2 --(p + a + X \P B + ob)) . (49) 



Taking into account Eq. (}32|) for H 2 and introducing the conformal time dr] = a l dt we 

get 

j- (a X ) = ±^y/(C - PoK + V* + Po )(«X) 4 , (50) 
which in principle can be integrated in terms of elliptic functions. 



5 Evolution of the scalar field 

In the previous section it was possible to find analytic solutions for the scale factor since 
the equations of motion (fT5Jl - (fT7j) combined to yield an equation for H decoupled from 
the scalar field $ and B-brane matter. In fact, it is also possible to derive an independent 
equation for \I/ if we retain only one type of source (cosmological constants or matter) 
similar on both branes (r s = T) and take the Universe to be flat (k — 0). Then, after 
introducing a new time variable JT] 

dp = h c dt, h c = H + ^j, J- = /' , (51) 
the equations (fT3 j) - (fT7j) together with (J2Uj) give a decoupled "master equation" for \E', 



(I - *)"^- - 3(2 - r) ^-j - 2((4 - 6*) - (4 - 3r)(l - 9)W(V) ) ^~ 

+12(2 - T)(l - - 8(4 - 3r)(l - y) 2 W(^) = 0, (52) 



where for the cosmological constants (r = 0) 

1 + (1-*)^ 
1 + (1 - m) 2 ^ 

and for matter (r > 0) 



= iTn , (53) 



1 + (1 - *) , . ,, 
W(9) = )p0 B K ' -* ' 3r • (54) 

i + (i-*) 2 

To solve this master equation analytically is still a formidable task, but in order to get a 
qualitative picture of the evolution of \& we can study Eq. (JH2*j) as dynamical system. 
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Cosmological constants 



Interval of \I/ 



Allowed region 



1) 


o 


>o, 


a + o B 


> o 


< x 


< 1 




|y| 


< 


\2xy/l - 


;r 


2) 


a 


>o, 


a + a B 


< o 


< x < 1 




cr 
(JB 


|y| 


> 


|2arVl - 


- x\ 














< X 


< 1 


M 


< 


|2aVl - 


x\ 


3) 


a 


<o, 


a + a B 


> o 


< x < 1 


-i 


cr 
(JB 


|y| 


< 


|2s>/l - 


- x\ 














< X 


< 1 


la/I 


> 


|2arVl - 


x\ 


4) 


a 


<o, 


a + a B 


< o 


0<x 


< 1 




M 


> 


|2xVl - 


x\ 



Table 1: Regions on the (x = ^,y = plane allowed for dynamics with cosmological 
constants. 

5.1 Cosmological constant dominated Universe 

By defining the variables x = ty, y = ^f' we can rewrite the equation (J52)l for T = as a 
dynamical system: 

x = y 

3y 3 ((4-6ar) -4(1 -x)W(x))y 2 A M 

V = —tt 1 r + - r r — -3y + Ax(l-x)W(x). (55) 

y 4x 2 (l-x) 4x(l-x) y V ) \ ! \ ! 

The physically allowed regions of the phase space (x, y) are constrained by the Friedmann 
equation (|15p. which in terms of the new variables turns out to be 

_ k 2 (cr + (1 - x) 2 cr°) 
K - 3x ,_ ~ 2 s • (56) 



4x' 2 (l-a) / 

If we let a and a B to take positive and negative values, there are four possibilities, sum- 
marized in Table ^ Note, Eq. (j^Bj) also assures that the time reparametrization (jB*T|) is 
acceptable — within the borders of the allowed phase space t-time and p-time always run 
in the same direction. 

In order to understand how to read the phase portraits correctly, it is instructive to 
see how the velocities \I/ get mapped to by 



^ dm dp , 



It turns out that the boundary 



B ± : y = ±2x^1 - x (51 



corresponds to \& = ±oo, while y = corresponds to ^ = 0. The outer reaches, y = ±oo, 
which are only relevant when a + (1 — x) 2 o B < 0, translate into 



y\y =±00 = ±kJ--x(1 - x)(a + (1 - x) 2 a B ), (59) 
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thus we get a minimal bound on the speed of the relative brane motion. At the borders 
x = and x = 1 we have \I> = and the same holds for the points (x = 0,y = 0) and 
(x = l,y = 0) as well. This tells us immediately that the brane collision = 0) can 
take place only extremely softly at vanishing relative speed — > 0~), and also when the 
branes approach infinite separation (\& = 1), they essentially stop moving with respect to 
each other (4f — > + ). At the crossing points 

0± = (x=i-fr, y=J3 ( 1 _^z) < cz) (60) 

the ratio % becomes undefined, meaning that with our redefinition of time (J5Tj) at these 
values of \P all velocities \I/ get mapped to a single respective velocity 
The system possesses the following fixed points: 

• Pi = (x = 0, y = 0), a saddle point with repulsive and attractive eigenvectors 
tangential to the upper and lower boundaries B ± , respectively; 

• P2 = (x = 1, y = 0), a (complex) spiralling attractor, but here the trajectories have 
to obey also the boundaries of the allowed region; 

and in the case a > 0, a + o B < there is also 

• P 3 = (x=l + ^-,y = 0),a saddle point (having de Sitter expansion for the scale 

factor EBBED- 

The limit of general relativity, = 0, \I/ = 1, necessarily holds for all trajectories drawn 
to the right boundary x — 1. The phase portraits in the four cases are depicted in Figure 
HJ let us consider these one by one. 

In the case 1) of o > 0, a + a B > (top left on Figure HJ) the allowed phase space is 
restricted into one bounded region, with \l/ = ±00 getting mapped to the boundaries B± 
and all the finite values being in the middle. Trajectories can either start near P 1; where 
the two branes very close to each other slowly start to move apart, or near P2, where 
the very far away branes slowly begin to move towards each other. The latter class of 
trajectories will eventually turn around, meaning these branes will never meet but travel 
back towards infinite separation. Hence all the trajectories are collected by the attractor 
P2, which corresponds to the limit of general relativity. (Of course, by definition the two 
fixed points provide two stationary trajectories themselves, i.e., two overlapping branes 
and two infinitely far away branes would always stay so.) 

In the case 2) with a > 0, a + a B < (top right on Figure ^) the allowed phase 
space consists of three areas, connected via two crossing points C±. Notice that in the 
left areas ^ = ±00 gets mapped to the boundaries B±, while y = ±00 corresponds to a 
finite velocity ([59)1 . Again, the trajectories can either start close to Pi or close to P2. The 
former trajectories explore the upper left allowed phase space region, pass through the 
upper crossing point C+ into the right allowed region and then proceed either to the fixed 
point P2, or pass though the lower crossing point CL into the lower left allowed region to 
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Figure 1: Phase portraits (x = \I/(p), y = ^'{p)) for cosmological constants on both 
branes: 1) a = l,cr s = -0.5 (top, left), 2) a = 0.5, o B = -1 (top, right), 3) a = 
-0.5, a B = 1 (bottom, left), 4) a = -l,a B = 0.5 (bottom, right). 

run into (x = 0,y = — oo), or in the marginal case terminate at the saddle point P3. These 
correspond to branes either starting very close to each other and moving apart towards 
infinity, or turning back and colliding, or, marginally, stopping at a finite distance. The 
other class of trajectories beginning near P 2 have the same possibilities, i.e., either they 
go back to P2, or pass through the crossing point C_ into the lower left region to run into 
(x — 0, y — —00), and there is also a marginal one that terminates at the saddle point P3. 
These correspond to branes starting to move closer to each other from a large distance, 
and then either separate again, or eventually collide, or in the marginal case come to a 
standstill at certain finite distance. 

In the case 3) with a < 0, a + a B > (bottom left on Figure HJ) the allowed regions 
are an inverse of the previous case. Now the common trajectories may only start near P2, 
explore the lower right region, necessarily pass through C_, visit the left allowed region, 
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pass through C + , and then either run to P 2 , or to (x — l,y — oo), or, marginally, to 
(x = l,y = 2/3). These correspond to far apart branes approaching each other to at least 
a certain distance, and then separating again towards infinity. All these trajectories will 
gently deliver the limit of general relativity (as they get to x = 1). The distinction where 
they end on the y axis is a coordinate artifact of time p, simply meaning how "fast" their 
speed measured in time t, i.e., \1/ tends to zero. 

In the case 4) with a < 0, a + o B < (bottom right on Figure the allowed part of 
the phase space consists of two separate regions. The trajectories can either begin near 
P2 and go to (x = 0,y = —00), or start near Pi and end up at P2, (x = l,y = 2/3), or 
{x = l,y = 00). The former trajectories describe branes coming closer from far apart to 
a collision, while the latter trajectories describe branes starting very close to each other 
and moving apart towards infinity, leading to the limit of general relativity. 

Some qualitative features of the radion dynamics in the cosmological constants case 
were inferred previously by Kanno, Soda, and Sasaki (201 by considering the Einstein 
frame radion potential, and by Webster and Davis [21] m the Jordan frame moduli space 
approximation. Both papers also point out that the A-brane Hubble parameter does 
not diverge at the event of brane collision, thus leading to the picture that branes could 
pass through each other with possible interesting cosmological consequences. Here we 
observe that the relative brane velocity in the cosmic time slows down and vanishes at 
the collision. However, this result must be taken with considerable caution, since one 
expects the effective action to receive corrections near the collision thus changing the 
behavior of the system 



5.2 Radiation dominated Universe 



For radiation on both branes (r = T B = 4/3) the dynamical system of Eq. (p)2*|) reduces 
to 



x' = y 



V 3 (4:-6x)y 2 
J 4x 2 (l-x) Ax{l-x) V ' 1 ' 

The Friedmann constraint is similar to the cosmological constants case, 



h2 = K 2 (p+(l-x)» = K 2 (p + (l-Vl/ )Vg) |^j-4 ^ ^ 

3X ~ 4z 2 (1-:e)) 3X I 1 ~ Ax 2 \l-x)) 



«0 



where in the last step we used Eqs. fITH)) . (f2l1|) to express p and p B . The allowed regions 
are listed in Table |2l including the exotic possibility of p B < 0. The discussion of how 
\1/ is mapped to ^' follows along the same lines as before yielding analogous results, only 
the outer reaches, y = ±00, relevant when p + (1 — \I/o) 2 Po < 0> translate into 




*\ y=±00 = ±m/--x(1 - x)(p + (1 - ^ ) 2 po s ) ( - ) • (63) 



a - - 2 
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Radiation 



Interval of \I/ 



Allowed region 



1) p B > 

la) p B <0,p + (l-*o) 2 Po >0 

2) p B <0,p + (l-^o) 2 Po <0 



< x < 1 1 2/ 1 < |2xvT^x 

0<x<l \y\<\2x^l^x 
< x < 1 \y\> \2xVl - x 



Dust Interval of \1/ Allowed region 

1) p B > < x < 1 1 2/ 1 < \2xVT^x\ 

la) p B < 0, p + (1 - ^ ) 3/2 Po > < x < 1 < |2aVT^x| 

< 0, po + (1 - ^o) 3/2 p^ < < x < 1 - fl . ^ ^ |y| > |2x^T^ 

1-^_<x<1 \y\<\2xVT^\ 



P B 



Table 2: Regions on the (x = = plane allowed for dynamics with radiation and 
dust (assuming p > 0). 

The limit of general relativity is still approached on the line x — 1. However, the main 
new feature is that the system (JUTj) does not have any fixed points, since the force term 
has cancelled. 

In the case 1) with p and p B both positive, and also in the more exotic case la) with 
p B < 0, but po > (1 — \&o) 2 Po the allowed phase space covers a single region (Figure 121 
top left). The trajectories can either start near (x — 0, y = 0) or (x = 1, y = 0), but all 
of them end somewhere on the y = line. This corresponds to close branes separating 
or far away branes approaching each other, but eventually all of them slowing down and 
stopping at some finite distance. The whole dynamics is ruled by friction, which is also 
in accord with our observation in Subsection 14. 21 that after a few e- folds of expansion the 
scalar field remains approximately constant. 

In the other exotic case 2) with p B < and p < — (1 — \l/o) 2 Po there are two allowed 
regions (Figure El top right). The trajectories beginning near (x = 0,y = 0) run through 
the upper region into (x = l,y = 0), meaning close branes separating towards infinity 
and reaching the limit of general relativity. The trajectories starting near (x = l,y = 0) 
run through the lower region into (x = 0,y = — oo), meaning far away branes coming to 
collide. 

5.3 Dust dominated Universe 

Taking both branes to be dominated by dust (r = r s = 1) the dynamical system is given 
by 

x' = y 
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Figure 2: Phase portraits (x = ^(p), y = *&'(p)) for radiation: 1) p = l,p s = 1, \I/o = 0.1 
(top, left), 2) p = 0.5, Pq = —l,^o — 0.1 (top, right), and for dust: la) po = 0.5, po = 
-1, = 0.5 (bottom, left), 2) p = 0.7, p£ = -1, tf = 0.1 (bottom, right). 



3,3 i (( 4 - fl )-(l- fWlW _ 5 + 



8x 2 (l-x) 4x(l-x) 2 

The allowed regions are determined by the Friedmann constraint, 

K 2 (p + (1 - X)V) _ « 2 (po + (1 - Xfl 2 {l - ^o) 3/2 Po S ) 



a 



a 



4x 2 (l-x) , 

where again in the last step we used Eqs. (fTHj). (|2*0j) to express p and p s . The regions 
available for the dynamics are listed in Table El Translating \1/ into proceeds with 
analogous results, only only the outer reaches, y = ±oo, relevant when p + (l — l I / o) 3 ^ 2 Po < 
0, correspond to 



V\ y=±00 = ±«y-|x(l - x) (po + (1 - x)V2(i _ ^ )3/2 p «) (^)~ 3/2 . (66) 
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The limit of general relativity is x = 1 as in the previous cases. The system (J64j) is 
endowed with three fixed points: 

• Pi = (x = 0, y = 0), a saddle point with repulsive and attractive eigenvectors 
tangential to the upper and lower boundaries B-t, respectively; 

• P 2 = (x = 1, y = 0), a (complex) spiralling attractor for p B > 0, but a saddle point 
for Pq < 0, notice that around here the trajectories have to obey also the boundaries 
of the allowed region; 

and in the case p$ < 0, p + Po (1 — \I/o) 3 ^ 2 > there exists also 

• p 3 = ( x = 1 - p ° 2(1 ~* o)3 , y = 0), an attractor. 

In the case 1) with p > 0, Pq > the phase portrait is qualitatively identical to case 
1) with cosmological constants (Figure^ top right). The trajectories can start either at 
Pi or P 2 , but they are all attracted to P2, approaching the limit of general relativity. The 
case la) with p B < 0, p + (1 — \I/ ) 3 / 2 Pq > offers novel features, since P 2 is now a saddle 
point, and the dynamics is ruled by the attractor P3 instead (Figure El bottom left). The 
branes do still start close by (Pi) or very far apart (P 2 ) but they all eventually move 
to a finite separation (P3), determined which does not bring about the limit of general 
relativity. In the case 2) with p B < 0, p + (1 — \E' ) 3 / 2 Pq < there is no attractive fixed 
point (Figure El bottom right). Close branes initially separate (Pi), but then meet again 
at a collision. 

Let us point out that taking the B-brane matter to be zero, i.e., considering a scalar- 
tensor theory with a single usual matter source, corresponds to the case 1) of radiation 
domination and case 1) of matter domination, yielding qualitatively identical phase por- 
traits in the respective cases. These can be compared with the investigations in general 
scalar-tensor cosmology ^Hl El E] carried out in the Einstein frame and using an inter- 
pretation of the scalar field in the "master equation" as a coordinate of a fictitious particle 
with velocity dependent mass. Our results agree - in the radiation domination regime the 
scalar field slows down and stops due to friction, while in the matter dominated regime it 
experiences attraction towards general relativity. However, the phase space method deliv- 
ers an easy to understand global picture without making any approximating assumptions 
and invoking the unphysical Einstein frame. 

6 Observational constraints 

6.1 The Eddington parameters and Solar system experiments 

The present state of the cosmological scalar field can be detected in the Solar system 
experiments [13j. Assuming that the B-brane matter has no significant influence on 
these experiments and taking into account Eq. (JHJ) for ufif), the observed values of the 
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gravitational constant G and the Eddington parameters 7, (3 read 

*l (4 + 2.„) 4^ 
*o(3 + 2u ) 3*o 

7 -l B -_ !_ = - 2 ' 1 -*'), (68) 

A 1- * 2 ^ggj _ 2(l-^o (69) 

^ 87rG(4 + 2cu )(3 + 2^o) 2 (4 - ^ ) 2 ' 1 ' 

Here \l/o denotes the cosmological value of the scalar field asymptotically far from Sun 
and co>o = u(^ ). 111 the case \l/o — 1 we get the general relativistic values (3 = 1 = 7 and 
87rG = k 2 . But in general $ 7^ 1 an d observational limits for 1 — ^ can be determined. 

The Lunar Laser Ranging experiment (LLR) [2H] gives the constraint 4/5 — 7 — 3 = 
(—0.7 ± 1.0) x 10~ 3 and the Shapiro time delay measured by tracking of the Cassini 
spacecraft [3U] gives 7 — 1 = (2.1 ± 2.3) x 10~ 5 . In our model 1 — \l/o > and this implies 
4/? — 7 — 3>0,7 — 1 < 0. It follows that in this case the observational bounds are in fact 
4/3-7-3 < 0.3 x 10~ 3 and 1 - 7 < 2 x 10~ 6 (cf. HE]). They yield observational limits for 
the asymptotic value of the scalar field and the modification of the gravitational constant 
correspondingly 

1 - *o < 10~ 4 (LLR), 1 - ^0 < 3 x 10~ 6 (Cassini) ; (70) 

k 2 /87tG < 1 - - x 10~ 4 (LLR), k 2 /8ttG < 1 - 4 x 10" 6 (Cassini) . (71) 
3 

The constraint on the present variation of is even more stringent. From the general 
expression [T3] 

^ = 3 + 2^q / 2a/Qgo) \ _ 4vj/ 

and observational data ^3] \G/G\ < 10 _11 yr _1 we get 

l^o/^ol < 10-V" 1 - (73) 

We can conclude that the asymptotic background for Solar system experiments is very 
near to general relativity — 1, ^ = 0. 



6.2 Nucleosynthesis constraints 

Primordial big-bang nucleosynthesis (BBN) provides a probe of a very early Universe and 
gives through primordial (relict) light element abundances (cosmologically relevant are D, 
3 He, 4 He, 7 Li) additional constraints for the parameters of the cosmological model under 
discussion. It can be used to test the non-standard cosmology as well as the non-standard 
particle physics. 

In the framework of general relativity and the standard model of particle physics 
the relict abundances depend on nuclear reaction rates Ri and only one cosmological 



19 



parameter, the ratio of baryon number density rij, to photon number density n 7 , 7710 = 
10 10 (^)- Here we will not discuss here the reaction rate uncertainties ARi, for a recent 
review, see |3~T] . Modifications of the standard cosmological model and the standard 
model of particle physics can be given in terms of two additional parameters: 1) the 
lepton asymmetry parameter (asymmetry between neutrinos and antinetrinos or neutrino 
degeneracy) £ e [H2]; 2) the speed-up factor = H(t)/ Hc(t) (in the framework of our 
model we considered it in Subsect. 4.2 for the radiation dominated Universe). 

Kneller and Steigman [33] have presented the primordial abundances {y D , Y p , y Li } of 
{D, 4 He, 7 Li} as linear functions of parameters {r?io, £, £ e } where £ = £,(t nuc i eosyn th)- The 
accuracy of these linear fits is approximately a few percent (lower than the current errors 
in observationally inferred value) over realistic ranges of parameters [321 EH] 4 < 7? 10 < 8 , 
0.85 < £ < 1.15 , —0.1 < £ e < 0.1 . The baryon to photon ratio r? 10 is determined at high 
precision from independent (non-BBN) measurements of CMB made by WMAP. In the 
standard ACDM theoretical model its value is 7710 = 6.14 ± 0.25 and it is approximately 
the same also in the case of nonstandard cosmological models |33]. We can take into 
account this (non-BBN) value and observed abundances of three elements, but then we 
cannot adjust the remaining two parameters {£, £ e } so that all three fits are satisfied. 
Since the lithium abundance derived from observations has possible problems, Kneller and 
Steigman (HSj derive possible values of the parameters only from D and 4 He abundances 
in two cases: 

• if £ e = 0, then rj 10 = 5.88t° f and £ = 0.969ig;gg| ; 

• if 7710 = 6.14 , then £ e « 0.037 and £ « 1.0203 . 

Now let us proceed with our specific scalar-tensor theory with coupling function (jHJ). 
The nucleosynthesis takes place at the beginning of the radiation dominated era and then 
the influence of the cosmological constant is negligible, so the relevant speed-up factor is 
given by Eq. (jUJ). The range of realistic values [33] 0.85 < £ < 1.15 translates for C/p 
as 

C 

- 0.28 < — < 0.32 (74) 
Po 

(small corrections from k 2 /8tiG ignored) and above-given observational constraints [HB] 
on £ yield 

• for £ ^ 0.969 we get C/p w -0.063 ; 

• for £ ^ 1.0203 we get C/p w 0.041 . 

The initial value of the dark energy density as a fraction of the background radiative 
matter energy density C/po is sometimes parametrized by the number of extra relativistic 
degrees of freedom, conventionally represented as additional neutrino flavors, AN U = 
(N v — 3) . In this case the speed-up factor reads 



f = J*L(l + LAN v ) — = —AN V . (75) 

5 8ttG V 43 J p A3 K J 
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Note that in our model the dark radiation is of geometric nature and its parametrization 
in terms of additional neutrino flavors is absolutely formal. 

Using updated predictions for light element abundances, BBN and CMB limit Cyburt 
et al. [HUIIIII] infer a combined limit —0.33 < AN U < 0.85, meaning 

C 

- 0.054 < — < 0.138. (76) 
Po 

This result can be compared with limits given by Ichiki et al. [36j for the ratio of dark 
and visible radiation in the framework of an one-brane model in an approximation p 2 « 0, 
since in this case if 2 and £ 2 coincide with the corresponding expressions in our two-brane 
model. Ichiki et al. ^ derived from BBN and CMB -0.41 < £ < 0.105 (using 
4.73 < 7/10 — 5.56 ). We see that limits for the positive dark energy are approximately the 
same, but the negative dark energy density is much more constrained by Eq. (J76|) . 

Finally let us note that the observational limits ([76^1 imply constraints on the param- 
eters of the cosmological model in the radiation dominated 

C = p + pg(l-fro) 2 _ i 

from Eq. ()45|) . For example, we can determine immediately that the case 2) with p B < 
0, po + (1 _ ^o) 2 Po — 0' discussed in Subsec. 15.21 is by far ruled out by the observations. 



7 Concluding remarks 

In this paper we have investigated a scalar-tensor cosmology arising as an effective 4- 
dimensional description of the 5-dimensional Randall- Sundrum two-brane cosmological 
scenario P in the Kanno-Soda low energy gradient expansion approximation The 
model is characterized by a scalar field $ (radion) measuring the proper distance of the 
branes and endowed with a specific coupling function (jHJ) implied by the construction of 
the model, a specific potential (fTTj) generated by cosmological constants on the branes, and 
barotropic fluid matter sources on both branes. In the limit when the (hidden) B-brane 
is empty, the model reduces to a scalar tensor cosmology with cosmological constant and 
matter, the limit of general relativity is obtained when \l/ = 1, \I/ = 0. We have assumed 
FLRW metric on both branes. 

Due to the specific form of the coupling function the Jordan frame equations of motion 
combine to give a scale factor dynamical equation decoupled from the scalar field and B- 
brane matter, leading to a Friedmann equation with only an additional dark radiation 
term reminding the extra dimension. This has enabled us to find exact analytic solutions 
for the flat (k = 0) Universe scale factor with various kinds of matter: cosmological 
constant, radiation, dust, and cosmological constant plus radiation. We have also found 
analytic solutions for the scalar field with radiation, and with cosmological constant plus 
radiation (latter case to be expressed in terms of elliptic functions). Over the years 
Barrow et al. have developed methods to obtain exact analytic solutions for scalar tensor 
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cosmologies with general coupling functions ^SlEHl- Forgetting the B-brane component, 
our model belongs to classes 1 and 3 studied in ref . [IS] . However, in the general case the 
solutions are not explicitly available in the cosmic time and one resorts to approximation 
techniques for early and late time behavior. Therefore, despite being a specific one, our 
model has an advantage of being directly solvable, thus providing an interesting example of 
a scalar-tensor cosmology throughout the entire radiation, dust and cosmological constant 
dominated eras. 

Although our model is governed by the Friedmann equation (with an extra dark ra- 
diation term) and hence in broad accordance with standard cosmology, we need to study 
the behavior of the scalar field, in order to be assured that the theory satisfies all obser- 
vational constraints and converges towards general relativity. A clever reparametrization 
of time ^1] allows to extract an independent evolution equation for the scalar field also 
when besides matter on the A-brane, there is similar matter (r s = T) on the B-brane. We 
examine the solutions of this "master equation" by applying the methods of dynamical 
systems. It seems determining the fixed points and drawing the phase portraits is a more 
straightforward approach than the usual procedure of changing into the unphysical Ein- 
stein frame and relying on an analogue with a fictitious particle with velocity dependent 
mass. Still, in the absence of B-brane matter, the results are in accord with previous 
general studies [101 E21 HI] : in a flat Universe with positive cosmological constant or dust 
the system relaxes to the limit of general relativity, while in the radiation case friction 
makes the scalar field to stop at an arbitrary position. 

Inclusion of the B-brane matter (possibly having negative energy density) complicates 
the story, as was also previously noted from the numerical studies in the moduli space 
approximation approach 8\. In summary, we find that general relativity, which in the 
branes' picture corresponds to infinitely far away branes, is the sole destiny for cosmo- 
logical constants a + a B > 0, or dust p > 0, p,f > 0- The solutions may approach 
general relativity or go the other way, i.e., to brane collision, depending on the initial 
conditions, for cosmological constants a + a B < 0, or radiation p + (1 — ^/ ) 2 p B < 0. 
Brane collision is the only future for dust with po + (1 — \&o) 3 ^ 2 Po — 0- Radiation with 
Po + (1 — \& ) 2 Po — has the radion stabilizing at arbitrary position due to friction. 
There are also fixed points in between the two extremes, namely a saddle for cosmological 
constants a > 0, o + a B < 0, and an attractor for dust p B < 0, po + (1 — ^o) 3 ^ 2 /^ > 0. 
The latter is especially interesting, since we see a natural stabilization of the radion for 
which one usually invokes an additional scalar field in the bulk [39 . 

The solar system tests restrict the present value of the radion to 1 — ^ < 10~ 4 , and 
allow it to be only very slowly varying, < 10 -11 yT -1 - The observations of D and 

He abundances together with CMB set limits on the fraction of dark radiation to normal 
radiation, —0.054 < ^ < 0.138. The results of these experiments can be accommodated 
in the model if its parameters are correspondingly constrained, which is not too hard, given 
that for large range of parameters general relativity is an attractor. Due to the specific 
coupling function (jHJ) the first integral (f2*3*j) for the two-brane cosmology coincides with 
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the Friedmann equation for the one-brane cosmology j2H] in the approximation where the 
matter density p on the brane is small and quadratic term p 2 can be neglected. This means 
that observational constraints for our model can be analyzed along the same patterns as 
in the case of the one-brane cosmology jSHl EZj- Also note that the specific coupling 
(|SJ) implies constant speed-up factor in the radiation dominated era and excludes a 
variable speed-up factor mechanism recently proposed by Larena et al. [IU] for solving 
the Li abundance problem. 
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